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Abstract 

In this paper we provide a formal derivation of both the Camassa-Holm equation and 
the fractional Camassa-Holm equation for the propagation of small-but-finite amplitude 
long waves in a nonlocally and nonlinearly elastic medium. We first show that the 
equation of motion for the nonlocally and nonlinearly elastic medium reduces to the 
improved Boussinesq equation for a particular choice of the kernel function appearing in 
the integral-type constitutive relation. We then derive the Camassa-Holm equation from 
the improved Boussinesq equation using an asymptotic expansion valid as nonlinearity 
and dispersion parameters tend to zero independently. Our approach follows mainly the 
standard techniques used widely in the literature to derive the Camassa-Holm equation 
for shallow water waves. The case where the Fourier transform of the kernel function has 
fractional powers is also considered and the fractional Camassa-Holm equation is derived 
using the asymptotic expansion technique. 

Keywords: Camassa-Holm equation, Fractional Camassa-Holm equation, Nonlocal 
elasticity, Improved Boussinesq equation, Asymptotic expansions. 


1. Introduction 


In the present paper we show that, in the long-wave limit, small-but-finite waves propa¬ 
gating in a one-dimensional medium made of nonlocally and nonlinearly elastic material 
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satisfy the Camassa-Holm (CH) equation [l| and the fractional CH equation (see equation 
(14. 101) ') when a proper balance between dispersion and nonlinearity exists. 

The CH equation 


Vr + KiV(^ + ivV(^ — = K2{2V(^V(^(^ + 


( 1 . 1 ) 


was derived for the propagation of unidirectional small-amplitude shallow-water waves 


ii,Qy,y,QQ 


when the nonlinear effects are stronger than the dispersive effects. Due to 
the fact that, even for smooth initial data, the solution of the CH equation stays bounded 
as its slope becomes unbounded, it is often used as an apprOTriate model capturing the 
essential features of wave-breaking of shallow water waves [7[. However, recalling that 
(HU is derived under the long wavelength assumption, it follows that the CH equation 
is valid only when the solutions and their derivatives remain bounded Q. For a discus¬ 
sion on a different criterion for wave-breaking in long wave models we refer the reader 
to It is interesting to note that the CH equation as a model for wave-breaking of 
water waves is an infinite-dimensional completely integrable Hamiltonian system [9|, [lOj . 
Another interesting property of the CH equation is the existence of the so-called peakon 


solitary wave solutions when ki = 0 [l|. At this point, it is worth pointing out that 
the derivation in the present study is also based on the long wavelength assumption 
and that ki is nonzero for the resulting equation. In addition to the studies about wa¬ 
ter waves, there are also studies that derive the CH equation as an appropriate model 


equation for nonlinear dispersive elastic waves. We refer the reader to 


ll| for the deriva¬ 


tion of a CH-type equation governing the propagation of long waves in a compressible 

hyperelastic rod, and to [ 1^ 1 for the derivation of a two-dimensional CH-type equation 

governing the propagation of long waves in a compressible hyperelastic plate. However, 

these studies relied only on the ’’geometrical” dispersion resulting from the existence of 

the boundaries, that is, from the existence of a bounded elastic solid, like a rod or a plate. 

Another type of dispersion for elastic waves is the ’’physical” dispersion produced by the 

internal structure of the medium. Therefore, one interesting question is to investigate 

whether the CH equation can be derived as an asymptotic approximation for physically 

dispersive nonlinear elastic waves in the absence of the geometrical dispersion. In this 

study, we consider the one-dimensional wave propagation in an infinite, nonlinearly and 

nonlocally elastic medium whose constitutive behavior is described by a convolution in- 
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tegral. We then show that, for an exponential-type kernel function, the CH equation 
can model the propagation of elastic waves even in the absence of the geometrical dis¬ 
persion. Furthermore, by considering a fractional-type kernel function we are able to 
derive a fractional-type CH equation, which indicates the possibility of obtaining more 
general evolution equations for suitable kernel functions. It is well-known that the KdV 
and BBM equations are valid at the same level of approximation while the CH equation 
is more accurate than the KdV and BBM equations. Therefore, when we neglect the 
highest order terms in the asymptotic expansion, the KdV and BBM equations and their 
fractional generalizations are also obtained as a by-product of the present derivation. We 
underline that the asymptotic derivation of the CH equation needs a double asymptotic 
expansion in two small parameters characterizing nonlinear and dispersive effects. How¬ 
ever, assuming simply that the two parameters are equal, the asymptotic derivations of 
the KdV and BBM equations can also be based on a single asymptotic expansion in one 
small parameter resulting from the balance of nonlinear and dispersive effects. 

The paper is organized as follows. Section 2 presents the governing equations of 
one-dimensional nonlocal nonlinear elasticity theory and gives the equation of motion 
in dimensionless quantities for various forms of the kernel function. In Section 3, using 
a multiple scale asymptotic expansion, the CH equation is derived from the improved 
Boussinesq (IBq) equation which is the equation of motion for the exponential kernel 
function. Section 4 presents the derivation of a fractional CH equation from the equation 
of motion corresponding to a fractional-type kernel function. 


2. A One-Dimensional Nonlinear Theory of Nonlocal Elasticity 


We consider a one-dimensional, infinite, homogeneous, elastic medium with a nonlin- 


[ 13 , 14, 


15l| and the references cited therein 


ear and nonlocal stress-strain relation (see 
for a more detailed discussion of the nonlocal model). In the absence of body forces the 
equation of motion is 

PqUu = {S{ux))x, (2.1) 

where the scalar function u{X,t) represents the displacement of a reference point X at 

time t, po is the mass density of the medium, S = S{ux) is the stress and the subscripts 

denote partial derivatives. In contrast with classical elasticity, we take the constitutive 
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equation for the stress S' as a general nonlinear and nonlocal function of the strain ux- 
That is, we assume that the stress at a reference point is a nonlinear function of the 
strain at all points in the body. As in [l^ the constitutive equation of the present 
model has the form 


S(A, t)= [ ai\X - Y\)a{Y, t)AY, a{X, t) = W'{ux{X, t)) 

Js. 


( 2 . 2 ) 


where a is the classical (local) stress, W is the strain-energy density function, Y denotes 
a generic point of the medium, a is a kernel function to be specified below, and the 
symbol ' denotes differentiation. The kernel a acts as a weight function that determines 
the relative contribution of the local stress cr(Y, t) at a point T in a neighborhood of X 
to the nonlocal stress S{X,t). So, when the kernel becomes the Dirac delta function, 
the classical constitutive relation of a hyperelastic material is recovered. Assuming the 
reference configuration is a stress-free undistorted state, we require that W{0) = W'(Q) = 
0. We point out that if we take W{ux) = (A -|- 2/r)(Mx)^/2 where A and fi are Lame 
constants, the above e quat ions reduce to those of the linear theory of one-dimensional 
nonlocal elasticity (see 13|). 

Without loss of generality, for convenience, the strain-energy density function may be 
considered to consist of a quadratic part {uxY/2 and a non-quadratic part G{ux) with 
G(0)=G"(0) = 0: 

W{ux) = l[]^{uxf + G('ux)], 

where 7 is a constant with the dimension of stress. Differentiating both sides of (EH) 
with respect to X and using (12.21) we obtain the equation of motion for the strain: 


PQUxtt = 7{ / ct{\X - Y\)[ux + g{ux)\AY}xx, 

Jr 


(2.3) 


where g{s) = G'{s) with g{0) = 0. Now we define the dimensionless independent variables 


A 


77 = - 

I V Po 


where I is a characteristic length and from now on, and for simplicity, we use u for ux 
and t for 77. Thus, E31) takes the form 


Utt = (/3 * (M-|- 5 (?x)))a:a;, 
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(2.4) 




where the convolution operator * is defined by 


P:^v= ^{x - y)v{y)dy 

Jr 


and P{x) = /a(|a;|). From a wave propagation point of view, the harmonic wave solutions 
to the linearized form of (12.41) are dispersive and the sole source of dispersion in the present 
model is the internal structure of the medium but not the existence of the boundaries. 
In general the kernel function 0 i s even, nonnegative and monotonically decreasing for 
a; > 0 (we refer the reader to [l3l| for the properties that an admissible kernel function 
must satisfy). A list of the most commonly used kernel functions is given in [l^. Here 
we consider two kernel functions: the exponential kernel [l^ which is the most widely 
used kernel function in the engineering applications of nonlocal elasticity QQ. and 
a fractional-type kernel function. These two kernels are chosen because they are the 
simplest representatives of the kernels that are convenient for asymptotic expansions. 
Moreover, as discussed in R,emark l4.21 starting from (IMI) with a general kernel satisfying 
some mild assumptions will lead to the same results with those of the two representative 
kernels. 

The exponential kernel is given by 0{x) = The Fourier transform of 0 is 

/3(0 = (1 + where ^ is the Fourier variable. Note that 0{x) is the Green’s function 
for the operator 1 — where represents the partial derivative with respect to x. 
Now, using the convolution theorem stating that the Fourier transform of the convolution 
of two functions is the product of their Fourier transforms, we take the Fourier transform 
of both sides of (12.41) . Then, substituting 0{^) of the exponential kernel into the resulting 
equation and taking the inverse Fourier transform, we obtain the equation of motion 
corresponding to the exponential kernel. Thus, for the exponential kernel, the equation 
of motion, (1231), reduces to the IBq equation 


'^xx '^xxtt — (^('^))a: 


(2.5) 


We next consider a fractional-type kernel function whose Fourier transform is 0{^) = 
(1 -I- (C^)"^)”^ where v may not be an integer. Note that the previous case corresponds to 
1 / = 1. To ensure the local well-posedness of the Cauchy initial-value problem defined for 


the resulting form of the equation of motion [l^, we impose the condition v >\. Noting 

that 0{x) is the Green’s function for the operator 1 -|- this time the equation of 
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motion, (EH), becomes an improved Boussinesq equation of fractional type 

“t" ( ^x) ~ {,9i.^y]xx ■ (^■^) 

Here the operator {—D"^Y is defined as {—DlYq = where T and 

denote the Fourier transform and its inverse, respectively. In Sections 3 and 4, starting 
from (ESI) and ( 1 ^ . respectively, and restricting our attention to quadratic nonlineari¬ 
ties, we will investigate the wave equations describing the unidirectional propagation of 
small-but-finite amplitude long waves. Knowing how the operator {—D'^Y scales under 
the scaling transformation X = dx where i5 is a positive constant will be important for 
the calculations in Section 4. To clarify the situation, let us consider the transformation 
q{x) = Q{X). Then the relationship between the Fourier transforms of q{x) and Q{X) 
is 

m) = [ e-^^^q{x)dx = i / e-^^'^Q{X)dX = iQ(|). (2.7) 

^ Jr 0 0 

Thus, using the inverse Fourier transform we get 

{-DlYq{x) = = 

= S^Y-DjcYQiX). ( 2 . 8 ) 

We also note that the operator {—D^Y is a translation invariant operator, so that 
{—D'^Yidix + 7 )) = {{—D‘Y)''q){x + 7 ) where 7 is a constant. 

3. Derivation of The Camassa-Holm Equation in the Long Wave Limit 

In this section we provide a formal derivation of the CH equation from (12.51) with g{u) = 
V? in the long wave limit. We restrict our attention to the quadratic nonlinearity since 
the KdV, BBM and CH equations all contain only quadratic nonlinearities. For a data 
on a compact support, (EH) has both right-going and left-going wave solutions that are 
moving apart. Assuming the two waves no longer overlap for a sufficiently large time, 
in the rest of this section we consider right- going, small-but-finite amplitude, long wave 
solutions of (12.51) . We first introduce the scaling transformation 

u{x^ t) = eU{6{x — t), 6t) = eU (Y, S) (3.1) 
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with Y = 6{x — t) and S = 6t to make the asymptotic behavior of (12.51) more transpar¬ 
ent. Here, small parameters e > 0 and (5 > 0 measure nonlinear and dispersive effects, 
respectively; more precisely, e denotes a typical (small) amplitude of waves whereas S 
denotes a typical (small) wavenumber, which is equivalent to supposing that the waves 
are long. Following the approaches in Q S, Q, Q we seek a solution in the form of a 
double asymptotic expansion in two small parameters e and S. The double limit process 
allows us to control dispersive and nonlinear effects in capturing the CH equation. Note 
that the asymptotic derivations of the KdV and BBM equations can be based on a single 
asymptotic expansion in one small parameter by taking, for instance, 5^ = e. Insert¬ 
ing the scale transformation (ED into (12.51) and multiplying the resulting equations by 
lead to the equation 


Uss - “^Uys - 6‘^{Uyyyy + Uyyss - ‘2Uyyys) - f-{U‘^)YY = 0 , 


(3.2) 


for U (Y, S). We now seek an asymptotic solution of (13.21) in the form 

C/(Y, e, S) = Uq{Y, S) + eUi{Y, S) + 5^2(Y, S) + S) + 0(e^ d^) (3.3) 

as e —>■ 0, (5 —>■ 0. Note that odd powers of 5 will never appear due to the existence of 
even order spatial derivatives only in the IBq equation. We assume that the unknowns 
Un {n = 0,1,2,...) and their derivatives decay to zero as |Y| —>• oo. Substituting the 
approximate solution (13.31) into (13.21) we obtain a hierarchy of equations in powers of e 
and i5^. The leading-order equation is given by 


{Ds - 2 Dy)Uos = 0 . 


(3.4) 


Recalling that we have already restricted our attention to the right-going solutions we 
take Uqs = 0 which implies Uq = Uo{Y). At 0(e), we have 


(Ds — 2 Dy)Uis — ((Uo)‘^)yy = 0 . 


(3.5) 


Differentiating this equation with respect to S we get an equation for Uiss which is the 
same as (13.4|) . With a similar argument, we take Uiss = 0. Then the solution of (13.51) is 
obtained as 

Uis = --^((Uo)‘^)y- 


( 3 . 6 ) 






At O(S^), we have 

{Ds — 2Dy)U 2S — UoYYYY = 0, (3.7) 

where we have used the fact Uqs = 0. Differentiating this equation with respect to S 
and then following the same line of reasoning, we take U 2 SS = 0- Then the solution of 
(13.71) is 

U 2 S = —-jUoYYY- (3.8) 

Finally, at C>(ei5^), the equation to be satisfied is of the form 

{Ds — 2 Dy)U3s — Uiyyyy + ^Uiyyys — 2.{UoU2)yy = 0, (3.9) 

where the fact Uiss = 0 has been used. Differentiating this equation twice with respect 
to S we get an equation for Uzsss which is the same as dSH). To follow the right-going 
wave, again we take U^sss = 0. Differentiating (EH) with respect to S and using dSH, 
(15^ and Ussss = 0, we get 

Usss = -^iiUo)^)YYYY + ■^{UoUoyyy)y- (3.10) 

Substituting this result into dSH gives Uss in the form 

Uss = —{UoU 2 )y — -^Uiyyy — ■^{UoUoy)yy + -^UoUoyyy- (3-11) 

Combining the above results with (13.31) we obtain 

Us = eUis + S^U 2 S + e 6 ^U 3 s + 0 {e^S^) 

= [{Uof + 2S^UoU2]y- y {Uo + eUi)yyy 

— -^{^UqyUqyy + 2.UqUqyyy) + 0{e^,S‘^). (3-12) 

At this level of approximation, adding Uqs = 0, (13.61) and (13.81) to (I3.12p . we form the 
equation for U{Y, S; e, <5) as follows 

Us + cUUy + —Uyyy H—y (SDyLVv -I- 2UUyyy) = 0. (3.13) 

We note that this equation is not in the standard CH form, for instance, the term Uyys 
of the CH equation is absent. But this can be remedied if we rewrite (13.131) in a moving 
frame defined by 


X = aY + hS, T = cS 
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(3.14) 






where a, b and c are positive constants. The free parameters a, b and c will be chosen later 
so that (|3.13l) becomes the CH equation in the new frame defined by X and T. Rewriting 
(13.131) in the new coordinate system with U{Y, S) = V(aY + bS, cS) = V(X, T) we obtain 


cVt + bVx + iCiWx + 




eS^a^ 


Vxxx + 


■(9VxVxx + ‘2VVxxx) = 0. 


(3.15) 


This equation still lacks a term corresponding to z;rCC dl.R . To remove the term Vxxx 
in favor of Vtxx , as in QflB , we use 


Vxxx = —jVtxx — -^^iyVx)xx ~\- ,e5^), (3.16) 

obtained from (I3.15L which is valid at this level of approximation. Thus, (13.151) takes 
the following form 

Vt+-^x +—VVx —^Vtjvx+ [3(3——)Vykxv:+2(l —-)RVxji:x] =0. (3.17) 

To remove the parameters e and 5 from this equation we introduce the scaling transfor¬ 
mation 

v = eV, X = 5C, T = 5t. (3.18) 


Then ()3.17p becomes 


b a or or 2a, a. , „ 

Vr + -VC + -vvt: - igOrCi + —[3(3 - —)vqVqq Y 2(1 - -)vvccc] = 0. (3.19) 

A special form of the CH equation given by dllD, which corresponds to the special values 
of Ki and K 2 , can now be recovered by making suitable choices for the parameters a, b and 
c. For a standard CH equation, the ratio of the coefficients of the terms and 


BBB (see B 


for a detailed discussion on 


is important and it must be 2:1, respectively 
this issue). In order to recover dnj we also impose the conditions that the coefficients 
of the terms vv,^ and Vr((; are 3 and -1, respectively. To ensure that all three conditions 
are satisfied we require 

•“Tr 

With this choice of the parameters, (13.171) and (13.191) reduce to 

6 

Vt + -^Vx + SeVVx — S'^Vtxx - -^{‘^VxVxx + VVxxx) = 0, (3-21) 
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and 


6 9 

Vt + “ '^CC'T = ^(2^^C'^CC + ^^^CCC)> (3.22) 

respectively. We note that (13.2211 is the standard CH equation (11.11) with m = 6/5 and 
K2 = 9/5. 

From (13.11) . (I3.14L (13.181) and (13.201) it follows that the coordinate transformation 
between (C,'r) and {x,t) is 


. 2,3, 2 

Q=—={x--t), T= — 

Vb 5 3^5 


(3.23) 


By applying this coordinate transformation, we can write (13.2211 in the original reference 
frame as follows 


3 5 3 

Xt Xx “b XVx "^Vxxx "^Xxxt — ~^\^XxVxx “b XVxxx) 


(3.24) 


with V = v{x, t). 

Remark 3.1. We point out that if we neglect the terms of order eS^ and higher in O 
we obtain the Benjamin-Bona-Mahony (BBM) equation ]1 a/ , instead of US. 2 ‘A) . In other 
words, if we seek an asymptotic solution of id. 2)) in the form 

U(Y, S] e, 6 ) = Uo{Y, S) + eUi{Y, S) + d'^U^iY, S) + 0(e^ ed^ <5^) (3.25) 

and if we follow similar steps in this section, we get the Korteweg-de Vries (KdV) 
equation 


Us + eUUy + -^Uyyy = 0 

instead of Sd.ltA) and then the BBM equation 

Vr + KlV(^ + 3vVi^ — = 0 , 


(3.26) 


(3.27) 


with Ki = Zaf j2. (where a is an arbitrary positive constant), instead of id.22]) . In the 
original reference frame (x,t), id. 2t)f) and ld.2''I\) are written in the form 


and 


Vt+Vx+ VVx + -Vxxx = 0 


3 5 „ 

H" H" . '^xxx . ‘^xxt — 0 

4 4 


(3.28) 

(3.29) 


with V = v{x,t), respectively. 
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4. Derivation of a Fractional Camassa-Holm Equation 


In this section we will handle (1^ with g{u) = and explore how to extend the 
asymptotic expansion of the previous section to (EH). As we did in the previous sec¬ 
tion, to make the asymptotic behavior of (12.61) more transparent, we use the scaling 
transformation EH in EH and we get 

Uss — ‘2 Uys + 5^'^{—DyY{Uss + Uyy — ‘^Usy) — = 0. (4.1) 

We then seek an asymptotic solution of (14.11) in the form 

U{Y, S- e, 5) = Uo{Y, S) + eC/i(y, S) + 6 ^’'U 2 {Y, S) + eS^’^U^Y, S) + 0(e^ <5^") (4.2) 

with the assumptions we made for EH- Substitution of EH into EH leads to a 
hierarchy of equations in powers of e and . As the process is quite similar to that of 
the previous section, we outline the basic steps involved in the asymptotic derivation. 
Again, a solution of the leading order equation is Uq = Uo{Y). Similarly, at 0(e), 
the asymptotic equation admits the solution given by (13.61) . At the asymptotic 

equation admits the solution 

U 2 S = li-D^YTUoY. (4.3) 

At 0(e(5^'^), the equation to be satisfied is of the form 

(Ds — 2 Dy)U3s + {—Dy)'' (UiYY — 2t/isv) — 2.{UqU2)yy = 0, (4.4) 

instead of EH- Following the same line of argument we obtain the solution of (EH as 
U^s = -{UoU2)y + \(-Dl)'^UiY + \{-Dl)YUe>UoY) - \uo{-D^y)’'UoY. (4.5) 

Substituting the above results for Uis, U 2 S and U^s into 

Us = eUis + Y’^U 2 S + + 0 {e^, <5^") (4.6) 

we get 

Us + eUUY - —{-DIYUy - —[?.{-DIY{UUy) - U{-D^YUy] = 0. (4.7) 
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for U{Y,S;€,S), which is valid at this level of approximation. Using the coordinate 
transformation (13.141) , we rewrite (14.71) in the new coordinate system with U(Y, S) = 
V{aY + bS, cS) = V{X, T) and then remove the term {—D\yVx in favor of {—D\YVt'- 

Ut + -Ux + -UUx + 1 i-DjcYVT 

c c Zb 

+ [(y - 2,){-Dlr{VVx) + V{-DlrVx\ = 0. (4.8) 

Using the scale transformation (13.181) and following the same line of argument that led 
to (13.211) and (|3.22l) . we obtain 

Ci O Q 

Ut + -Ux + -^{VYx + S^Y-Dj.rVr + -^[2{-D%y {VVx) + U(-7?|)"Ux] = 0, 

(4.9) 

instead of ()3.2ip and the fractional CH equation 

Vt + + i-D^Yvr = -^[2{-Dly{vvy + vi-D^Yvi^], (4.10) 

instead of (13.221) . In the present case, the parameters a, b and c are obtained as 

a=(^)i/2^^, b=^a, c=^. (4.11) 

To the best of our knowledge, the fractional CH equation, (I4.10|) . has never appeared 
before in the literature. 

From (13.11) . (I3.14|) . (13.181) and (14.111) we get the coordinate transformation between 
(C, t) and {x, t) as 

C = a{x-\t), T=^t. (4.12) 

5 3 

Using this coordinate transformation, (jd.lOp can be rewritten as 
Vt+v^+vv^ + ‘^{-DlYv^ + ^{-DlYvt = -^[2{-DlY{vvY} Y v{-DlY (4.13) 
for V = v{x,t) in the original reference frame. 

Remark 4.1. As in the previous section, if we neglect the terms of order eS^'^ and higher 
in we obtain a fractional BBM equation, instead of That is, if we seek an 

asymptotic solution of in the form 

UiY, S- e, 5) = Uo{Y, S) + eUi{Y, S) + 5‘^''U2{Y, S) + 0(eYeS^Y 
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(4.14) 





























and if we follow similar steps in this section, we get the fractional KdV equation ]2(\I 


Us + eUUy - ^—{-DIYUy = 0 


instead of o and then the fractional BBM equation \2m 

Vr + + ivVQ + { — D^YVr = 0 , 


(4.15) 


(4.16) 


with Ki = j2 (where a is an arbitrary positive constant), instead of |^.10[ ). In the 
original reference frame (x,t), {4-15^ and take the following forms 

1 


and 


Vt+Vx+ VVx - -^i-D(,Yvx = 0 


Vt+Vj:+ VV:c + ^ (--D^)''^'a; + ^{-D^Yut = 0 


(4.17) 


(4.18) 


with V = v{x,t), respectively. 


Remark 4.2. In Sections 3 and 4, we have derived the CH, KdV and BBM equations and 
their fractional counterparts by starting with g and considering the kernels defined by 
and P{f,) = (l + {((^Y) ^ ■ ^^6 natural question is what would happen 
if we had started with an arbitrary kernel function f3{x). To provide the answer to this 
question, let us suppose that the Fourier transform 13(f) is regular enough to have an 
expansion of the form, up to constants, 

i/m = ^ + {eY+o{ie)^) 


with 1 < V < p for f near 0. When we consider an asymptotic expansion (like 
for a solution of we observe that the term O ((C^)^) in the above equation yields 

to O (Y^) terms. Thus, only the leading-order terms 1 + (Y)'^, corresponding to the 
operator 1 + {—DY)^, will be effective at the level of approximation considered in this 
study. In other words, up to constants, the asymptotic derivation will give rise to exactly 
the same set of evolution equations, that is, i4-35\) and (|4.16‘| ). In that sense our 

results in Sections 3 and 4 in fact hold for with a more general class of kernels. 
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